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ABSTRACT. Remarkable parallelism between the theory of integrable systems of 
first-order quasilinear PDE and some old results in projective and affine differen- 
tial geometry of conjugate nets, Laplace equations, their Bianchi-Backlund trans- 
formations is exposed. These results were recently applied by I.M.Krichever and 
B.A.Dubrovin to prove integrability of some models in topological field theories. 
Within the geometric framework we derive some new integrable (in a sense to be 
discussed) generalizations describing A-wave resonant interactions. 



Ten years ago [10] a natural hamiltonian formalism was proposed for the class of 
homogeneous systems of PDE 
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(called "one-dimensional systems of hydrodynamic type"). Later (see [11]) it was 
generalized for the class of multidimensional 
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u % = u l (x, t), i = 1, . . . , n 
systems. 

Some systems (1+1, h) of physical importance such as Whitham equations (the 
averaged 1-phase KdV equation) and Benney equations have the notable property of 
being diagonalizable: under a suitable choice of field variables u % (Riemann invariants) 
the equations become 

(1) ul(x) = Vi(u)ui 

(there is no summation over As we have proved in [32], these properties (hamilto- 
nian property and diagonalizability) imply integrability. Deeper insight into this type 
of integrability is given by the theory of orthogonal curvilinear coordinate systems. 
This chapter of classical differential geometry was being intensively developed at the 
beginning of the XX century ([6], [8], [20]). In fact this theory gives the geometric 
background for integrability of systems (l+l,h), (N+l,h), (1+1, nh). These forgotten 
corners of differential geometry seem to be worth revisiting. 

An example. The well-known Bullough-Dodd-Jiber-Shabat equation u xx — u u = 
e u — e~ 2u (in the form (In h) uv = h — l/h 2 ) was introduced for the first time in [35] 
where the respective linear problem was given as well as a proper Backlund trans- 
formation for it! It is much simpler and "geometric" than Backlund transformations 
discussed recently [1], [31] in the context of integrable systems. 

In this paper we will sketch some applications of methods originating from classical 
differential geometry to equations of types (1+1, h), (N+l,h), (1+1, nh). 

1. Diagonal systems of hydrodynamic type and orthogonal curvilinear 
coordinate systems in R n 



Let us recall briefly the main results of [10], [32]. A (generally nondiagonal) system 
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v l j(u)u 3 x is hamiltonian if there exist a hamiltonian H — J h(u) dx and a 
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hamiltonian operator 

& = gV{u)j- + W{u)u k x 



which define a skew-symmetric Poisson bracket on functional 
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satisfying the Jacobi identity and generate the system 

PjTT 

(2) ul(x) = {«*(*), = = (9 ij d k djh + b*d s h)u k x = vi(u)u k x 

where d s = d/du s . B.A.Dubrovin and S.P.Novikov [10] proved that the necessary and 
sufficient conditions for A i3 - to be a hamiltonian operator in the case of non- degeneracy 
of the matrix g %d are: 

a) g l i = gi % , i.e. the inverse matrix g" 1 defines a Riemannian metric. 

b) b k J = —g ls Y 3 sk for the standard Christoffel symbols T 3 sk generated by 
9ij • 

c) the metric has identically vanishing curvature tensor. 

In such case we have Vj(u) = V l V jh = g ls V S V jh with the covariant derivatives 
defined by gij. 

Lemma [32], [33]. In order that a matrix Vj(u) be a matrix of a hamiltonian 

system (1+1, h) with a nondegenerate metric in A 13 it is necessary and sufficient that 
there exists a nondegenerate zero curvature metric g^ such that 

a) gikVj = g jk Vi and 

b) V jV l k = Vfci>], where V is the covariant differentiation generated by the metric 

gij- 

For a diagonal matrix Vj(u) = Vj(u)5 l j this implies that (see [32], [33]) g^ is also 
diagonal and 

(3) , 9iVk =r* =kln gifcfc , d i = d/du i 

[Vi - v k ) 2 

(hereafter we do not imply the summation on repeated indices!). From (3) we deduce 

(4) = a i+i + k . 

From a differential geometric point of view, to give a zero curvature nondegenerate 
diagonal metric is equivalent to giving an orthogonal curvilinear coordinate system on 
a flat (possibly pseudo- Euclidean) space (see [6]). Locally these coordinate systems 
are determined by n(n — l)/2 functions of two variables (L.Bianchi). A striking fact 
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can be discovered: formula (3) was found in [6] (p. 353)! This formula is crucial for 
the integrability property of diagonal hamiltonian systems (1): if we interpret it as 
an overdetermined (compatible in view of zero curvature property of g) system on n 
unknown functions Vj(u) (ga given) we can generate from every its solution Vj(u) a 
symmetry (commuting flow) 

u\ =Vi(u)u l x , i = l,...,n, 

of (1) and a solution of (1) (the generalized hodograph method, see [33] for the 
details). One can prove ([33]) the completeness property for this class of symmetries 
and solutions parameterized by n functions of 1 variable - the generic Cauchy data 
for our diagonal system (1). 

The corresponding geometric notion used in the theory of orthogonal curvilinear 
coordinate systems corresponding to (3) is the so called Combescure transformation 
(see [6]). 

Definition . Two orthogonal curvilinear coordinate systems x % = x l (u 1 , . . . , u n ) 
and x l = x'^u 1 , . . . , u n ) in the same flat (pseudo) Euclidean space R n = {(a; 1 , . . . , x n )} 
are said to be related by a Combescure transformation (or simply parallel) iff their 
tangent frames e*j = dx/du % and e*j = dx/du % are parallel in points corresponding to 
the same values of curvilinear coordinates u l . 

Let us take the quantities Hi(u) = |ei| = y/gTu Hi(u) = |efj| (Lame coefficients). 

Proposition The quantities Vi(u) = Hi(u)/Hi(u) satisfy (3) with = diH k /H k , 
the connection coefficients for the metric ga = Hf. Conversely for any solution Vi of 
(3) ga = (viHi) 2 will give an orthogonal curvilinear coordinate system related to the 
coordinate system with the metric ga = Hf by a Combescure transformation. 

The theory of Combescure transformations coincides with the theory of integrable 
diagonal systems of hydrodynamic type. 

Physical examples of such systems ( Whitham equations, Benney equations) have 
hamiltonian structures (2) with diagonal metrics ga possessing the so called Egorov 
property: dig kk = d k gu- As we have demonstrated earlier ([34]) this is a consequence 
of Galilei invariance of the original systems. See also [12] for the algebro-geometric 
background of this property for averaged integrable systems. Using this property and 
homogeneity of coefficients one can find explicit formulas for solutions of (3) for the 
systems in question [33], [34]. 

The class of Egorov orthogonal curvilinear coordinate systems is interesting in 
itself and merits our special attention. 

2. Egorov coordinate systems, the iV-wave problem and its generalizations 



Introducing j3i k (u) = diH k /Hi, i ^ k, (3u(u) = (rotation coefficients of a given 
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orthogonal curvilinear coordinate system with gu = H 2 , see [6]) one can easily check 
the following: 

a) vanishing of the curvature tensor is equivalent to 

(5) djfiik = PijPjk, j ^ k, 

(6) dipik + d k p ki + PsiPsk = 0, i^k. 

b) the Egorov property d; t gkk = dk9a reduces to 

(7) Pik = P ki- 
ln the Egorov case condition (6) is equivalent to TPik — 0, T — d\ + . . . + d n . Con- 
sequently the problem of classification of Egorov coordinate systems is reduced to 
description of all off-diagonal symmetric matrices (Pik) satisfying (5) and Tp ik = 0. 

B.A.Dubrovin [12] have recently observed that this problem coincides with the 
purely imaginary reduction of the well-known integrable system describing resonant 
TV- wave interactions. Namely, restriction of Pik on any (x, t) plane u l = a l x + b l t gives 
(compare, for example, [28]) 

[A,T t ]-[B,T x ] = [[A,T],[B,T]], 

A = diag^a 1 , . . . , a n ), B = diagib 1 , . . . , b n ), T = (p ik ) with additional reduction 
Im T = 0, T T = T. For the case N = 3 this reduces to 



(8) 



' b] + C x b\ = Kb 2 b 3 , 

b 2 t + c 2 b 2 x = nb 1 ^, 
_ b 3 t + c 3 b 3 x = Kb x b 2 . 



This is a system of type (1+1, nh), integrable by the 1ST method ([28]). 

Now we can compare the progress achieved in the modern integrability theory 
for (8) and the results obtained more than 70 years ago in the theory of Egorov 
coordinate systems initiated by G.Darboux in 1866 and continued by D.Th. Egorov in 
1901 in his thesis (see [13]). It was Darboux [6] who proposed to call this special type 
of coordinate systems Egorov type systems. From the point of view of integrability 
properties remarkable progress was achieved by L.Bianchi in 1915 (see [2]). He found a 
Backlund transformation for this problem and established the permutability property 
as well as the superposition formula for it. We shall remark here that the pioneering 
results on Backlund transformations and their permutability in the well-known theory 
of constant curvature surfaces in R 3 are due to Bianchi also. 

Let us take an orthogonal curvilinear coordinate system (not necessary of Egorov 
type) with Lame coefficients H^u) and rotation coefficients Pikiu). Bianchi applied 
to it a generalization of Ribaucour transformations known in the theory of transfor- 
mations of surfaces. We recall that two surfaces x(u, v) and x(u, v) in R 3 are related 
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by a Ribaucour transformation iff there exists a two-parametric family of spheres 
S(u, v) such that each sphere S(uq, v q) is tangent to both surfaces in corresponding 
points x(u ,Vo),x(u ,v ) and this correspondence preserves the curvature lines on 
the surfaces. For the case of a pair of orthogonal curvilinear coordinate systems in 
R 3 we need a three-parametric family of spheres (or an n-parametric family for the 
n-dimensional case) tangent in the corresponding points to one of three families of co- 
ordinate surfaces as well as to a coordinate surface of the other curvilinear coordinate 
system. Since due to the classical Dupin theorem coordinate lines in any orthogonal 
curvilinear coordinate system are curvature lines their correspondence is guaranteed. 
In terms of the rotation coefficients fii k one shall find a solution ji(u) of 

(9) da k = Pttfi, i ^ k, 

to define the corresponding Ribaucour transformation ([2]) 

(10) = p ik - ^(d klk + ]T (3 skls ), A = £( 7p ) 2 . 

71 s+k V 

For the case of Egorov systems we shall complete (9) and restrict 7i to satisfy 
Qili = C7i - S PiU" ( c = const ) 

or 

(11) T 7i = c 7i , f = d 1 + ... + d n 
Then the Backlund transformation in question is 

(12) hk = Pik - 

The permutability property for any orthogonal coordinate system requires a quadra- 
ture, but for Egorov systems it may be found explicitly and provides the following 
formulas for the fourth Egorov system /3i k related to obtained from a given 

Egorov system (3 ik with constants c', c" {d + c" ^ 0) and potentials 7^, r y'/ in (10): 



7i Ti / _; 1 _// \ / - ,/ \ 9 ' 1; 



One can enjoy reading [4], [23] where these formulas were rediscovered in the context 
of 3-wave system. So the basic integrability results for (8) were established long ago by 
Darboux, Egorov and Bianchi certainly with the exception of the 1ST transformation. 

An unexpected result (hidden in [6]) consists in existence of a homogeneous system 
(1+1, h) of three equations related to (8) by a nonlocal transformation. Geometrically 
this is trivial. Given an orthogonal curvilinear coordinate system in R 3 we have in 
each its point P(x ,y , z ) the orthogonal 3-frame of tangent planes 

(13) z = p k (x - xq) + q k (y - y ) + z , k = 1,2,3. 
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Let us parameterize it by 3 functions A(x, y, z), B(x, y, z), C(x, y, z), coefficients pk, qk 
of tangent planes being three solutions of 

J pq + Ap + Bq = 0, 

\ p 2 - q 2 + 2(Cp + Hq) = 0, 2{BC - AH) + 1 = 0, 

different from the trivial solution p = q = ([3]). Then the Frobenius compati- 
bility conditions for these three families of distributions (13) give a system of three 
homogeneous first-order equations of type (2+l,h): 



(14) 



' 2{AC Z - CA Z ) = 2C y + B y -A x 
2(BH Z - HB Z ) = 2H X + B y - A x 
AH Z - HA Z + BC Z - CB Z = A y -B. 



where H = (2BC - 1)/2A . 

For this system one can reformulate the Backlund-like transformation (10) given in 
terms of /3ik(u). A number of different transformations producing (with quadratures) 
solutions of (14) parameterized by arbitrary many functions of one variable may be 
found in [6]. Thus (14) is integrable in a sense to be discussed elsewhere. 

If we will search for solutions of (14) which do not depend on z then a remarkable 
integrable (1+1, h) system of three equation appears. Since one can easily prove the 
equivalence of ^-independence in (14) and the Egorov property (7) we have received a 
homogeneous system related to (8) by a nonlocal change of variables. In Euler tp, if), 9 
parameterization of orthogonal 3-frames it reads 

(ipt \ ( — cos 2 — sin (p cos <p/ sin 6 

9 t = — sin 9 sin (p cos (p — sin 2 (p 

tp t J \ — cos 0(1 + cos 2 tp) — sin (/? cos y?cos6>/sin# 1 

This nonlocal change does not affect the existence of higher order conserved densities. 
Recently Ferapontov [17] have proved the uniqueness result for such 3x3 homoge- 
neous systems possessing higher-order conserved densities: they may be transformed 
to (15) by reciprocal and point transformations. Also another nonlocal transition 
from (8) to (15) as given there. 

The matrix of (15) has constant eigenvalues — 1,0, +1 but its eigenvector fields 
(properly normalized) form so(3) Lie algebra, consequently (15) is a non-diagonalizable 
(l+l,h) integrable system. 

The complete system (14) certainly may be called a (2+l)-dimensional generaliza- 
tion of the (l+l)-dimensional 3-wave system (8). Orthogonal curvilinear coordinate 
systems in R n provide also only a (2+1) -dimensional generalization of the (1+1)- 
dimensional iV-wave system since they are parameterized by n(n — l)/2 functions of 
two variables (L.Bianchi). 




3. Semihamiltonian diagonal systems and coordinate systems with conju- 
gate lines 
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The class of integrable diagonal systems (1) is wider than the class of hamiltonian 
systems of this type. Namely, the property (4) which is a weaker consequence of 
the hamiltonian property is sufficient ([33]). Let us call a diagonal system semi- 
hamiltonian if n = 2 or if n > 2 and Vi(u) satisfy (4). As a physical example of a 
semihamiltonian (but non-hamiltonian for n > 3) system one can mention the ideal 
Langmuir chromatography and electrophoresis systems ([33]). 

To every semihamiltonian system we can relate a diagonal metric ga via 

di In g kk /2 = diV k / (vi - v k ) 

. This metric is not flat in general though some coefficients of the curvature tensor 
vanish as a consequence of (4). Namely introducing Hi = ^fgTi , Pi k (u) = diH k /Hi, 
we can find out that (4) is equivalent to the set (5) of equations on j3i k . Solutions of 
(5) may be parameterized by n(n — 1) functions of 2 variables. This system coincides 
with the compatibility conditions for a linear system 

dii) k = Pikipi, i^k. 

Restricting (5) on 3-dimensional planes u l = a l x + b' l y + c l z in R n we obtain (for 
general nonvanishing constants a\b\c l ) a (2+1, nh) system on n{n — 1) quantities 

(3 ik (x,y,z). 

As we have seen earlier the theory of hamiltonian diagonal systems (1) is closely 
related to the theory of orthogonal curvilinear coordinate systems in R n . The geo- 
metric background for the theory of semihamiltonian systems is given by the theory 
of coordinate systems with conjugate coordinate lines (see [6] and [7], t. 4, ch. 12). 
A general (non-orthogonal) coordinate system x(ui, u 2 , u 3 ) in R 3 is called a system 
with conjugate coordinate lines (or simply a conjugate coordinate system) if on ev- 
ery coordinate surface S io = {u % = const} in every point P(x , y , z ) the lines of 
intersection of this surface with two other coordinate surfaces belonging to other one- 
parametric families of coordinate surfaces and containing P(x , y , z ) are conjugate 
on Si (with respect to its second fundamental form). Every orthogonal curvilinear 
coordinate system is conjugate due to Dupin theorem mentioned above. The theory 
of conjugate coordinate systems was developed by Darboux and others and borrowed 
a lot of results from the classical theory of conjugate coordinate nets on surfaces in R 3 
(known as "nets" or "reseaux", see [14], [21], [22], [36]). A number of Backlund-like 
transformations for these coordinate systems was given with permutability properties 
(though the superposition formulas therein require quadratures). 

Every conjugate coordinate system x l = x^u 1 , . . . ,u n ) in R n = {(x 1 , . . . ,x n ) is 
characterized by the conditions of conjugacy of coordinate lines: 

(16) d l d k x = TUu)d k x + T\ k (u)d l x, i^k. 

This system of equations coincides with the system describing hydrodynamic type 
conserved quantities of a semihamiltonian system with ([33]). Quantities T ki in (16) 
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satisfy its compatibility conditions 

fy^ki = di^kji dj^ki = ^kj^ji + ^ki^ij ~ ^ki^kj^ '^J^' 

equivalent to the semihamiltonian property (4). Introducing Hi(u) as solutions of 
diHkiu) = r^Hk(u) and (3^ = diHk/Hi, i ^ k, we receive a set of satisfying 
(5). The converse is also true: given a solution fak of (5) one can find (a number 
of) semihamiltonian systems related to it. Any semihamiltonian system also may be 
related to a Combescure transformation of conjugate coordinate systems ([6]). 

This geometric interpretation provides another example of integrable (2+1, h) sys- 
tem. Namely, given a conjugate coordinate system in R 3 one can take the field of 
its (non-orthogonal) tangent 3-frames (e\, 62,63) and parameterize it by 6 indepen- 
dent functions e l k (x,y,z),i = 1,2, k = 1,2,3 (the coefficients ef, may be set to 1 
due to normalization). Then the Frobenius compatibility conditions give 3 homo- 
geneous first-order PDE on e' l k (x, y, z), i = 1,2, k = 1,2,3. Another 3 equations are 
given by the conjugacy condition det((e*j • V)efc,ej,efc) = 0,i < k. This system of 
6 equations is a homogeneous (2+1, h) system in question. Its ^-independent solu- 
tions satisfy a (l+l,h) system enjoying properties analogous to those of (15): it has 
constant eigenvalues — 1,0, +1 (all doubly degenerate) and 6 linearly independent 
fields of eigenvectors forming (if properly normalized) a nontrivial Lie algebra. This 
remarkable system will be studied in subsequent publications. 

4. Additional topics 

Recently O.I.Mokhov and E.V.Ferapontov [27], [16] found a nonlocal generalization of 
the hamiltonian formalism of hydro dynamic type (2). V.E.Ferapontov communicated 
to the author about the further generalization resulting in the following beautiful 
theorem: any semihamiltonian system (1) has a nonlocal hamiltonian structure with 
a hydrodynamic hamiltonian and a hamiltonian operator with (possibly infinitely 
many) nonlocal terms similar to those in [16]. Grinevich [19] derived a series of 
nonlocal symmetries for Whitham equations as well as original KdV equations. 

Weakly nonlinear semihamiltonian systems (i.e. systems (6) with d(Vi = without 
summation on i, such systems are also called "linearly degenerate") were studied in 
[15], [30]. The theory of such systems is connected to the theory of n-webs on Euclidean 
plane, Dupin eyelids and Stackel metrics (E.V.Ferapontov). Among the results are: 
quasiperiodic behavior of their solutions ([30]), complete description of such systems 
and complete sets of their hydrodynamic symmetries ([15]). E.V.Ferapontov com- 
municated to the author the following fact: any n-phase (n-zone) quasiperiodic (or 
a n-soliton) solution of the KdV equation can be represented with a solution of a 
weakly nonlinear semihamiltonian system R\ = (J2k^i R k )R % x -> i = 1? ■ ■ ■ ,n. These 
results shall be compared with Curro and Fusco's results [5] in the soliton-like inter- 
actions of Riemann simple waves for some 2x2 systems. 
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IST-like methods were developed in [18], [24], [25] for some diagonal hamiltonian 
systems of physical importance. Certainly this approach shall be related to our 
geometric methods. 

In a series of preprints (see [9], [26] and references therein) I.M.Krichever and 
B.A.Dubrovin exposed a remarkable link between the theory of Egorov coordinate 
systems and Witten-Dijgraagh-Verlinder-Verlinder equations for the correlation func- 
tions of topological conformal field theories proving their integrability. 

Acknowledgement. The author enjoys the occasion to thank Dr. E.V.Ferapontov 
and Professor J. Gibbons for many useful discussions and hospitality in Imperial Col- 
lege, London, where a part of this paper was written as well as Professor P. Clarkson 
and our colleagues in Exeter for their hospitality during the Workshop. 
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